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ABSTRACT 

We report the detection of a rare transiting brown dwarf with a mass of 59 Mj^p and radius of 1.1 Rjup around the metal-rich, [Fe/H] = -rO.44, G9V 
star CoRoT-33. The orbit is eccentric (e = 0.07) with a period of 5.82 d. The companion, CoRoT-33b, is thus a new member in the so-called brown 
dwarf desert. The orbital period is within 3% to a 3:2 resonance with the rotational period of the star. CoRoT-33b may be an important test case 
for tidal evolution studies. The true frequency of brown dwarfs close to their host stars (P < 10 d) is estimated to be approximately 0.2% which is 
about six times smaller than the frequency of hot Jupiters in the same period range. We suspect that the frequency of brown dwarfs declines faster 
with decreasing period than that of giant planets. 
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1. Introduction 

As of June 2015, W. R. Johnston lists 2085 confirmed and 562 
candidate brown dwarf^], 427 of them are marked as located in 
either a binary or multiple system. Only 65 brown dwarfs are 
companions to FGK dwarfs at orbital distances less than 2 AU 
from the star (Ma & Ge 2014). Of these only nine (and one addi¬ 
tional reported here) are transiting (Table 1) which provides the 
possibility of measuring directly the mass and radius directly. 

* The CoRoT space mission, launched on December 27th 2006, 
has been developed and is operated by CNES, with the contribution 
of Austria, Belgium, Brazil, ESA (RSSD and Science Programme), 
Germany and Spain. Based on observations made with HARPS (High 
Accuracy Radial velocity Planet Searcher) spectrograph on the 3.6- 
m European Organisation for Astronomical Research in the Southern 
Hemisphere telescope at La Silla Observatory, Chile (ESO program 
188.C-0779). 

** Based on observations obtained with the Nordic Optical Telescope, 
operated on the island of La Palma jointly by Denmark, Finland, 
Iceland, Norway, and Sweden, in the Spanish Observatorio del Roque 
de los Muchachos of the Instituto de Astrofisica de Canarias, in time 
allocated by the Spanish Time Allocation Committee (CAT). 

* See at http://www.johnstonsarchive.net/astro/browndwarflist.html I. 

The DwarfArchive.org project, maintained by C. Gelino, 
D. Kirkpatrick, M. Cushing, D. Kinder, A. Burgasser 

(http://spider.ipac.caltech.edu/staff/davy/ARCHtVE/index.shtml I 
tabulates 1281 L, T and Y dwarfs but it does not contain information 
about the masses and radii of the catalogued objects. The Johnston’s 
catalogue reports radius values for 191 brown dwarfs, but most of 
these were obtained by several different and often indirect methods, 
e.g. model-fitting to the spectral energy distributions. Radii of brown 
dwarfs measured by direct methods are listed in Table 1. 


Given the meager observational data, an increase in the number 
of brown dwarfs with accurate parameters is necessary to obtain 
a better understanding of the formation, structure, and evolution 
of these objects. Transiting brown dwarfs offer the best possibil¬ 
ity for characterizing them in a model-independent way. 

The radii of brown dwarfs and giant planets are in the same 
range around 1 Rjup, but the masses are quite different. Brown 
dwarfs straddle a border at 65 Jupiter masses. Objects with 
masses below this border fuse deuterium (D) while those above 
this mass limit fuse lithium (Li) in episodic events (e.g. Chabrier 
et al. 1996), and planets do not ignite their material. Spiegel et 
al. (2011) showed that the lower mass limit for brown dwarfs - 
defined by deuterium-ignition - is between 11-16 Mjup depend¬ 
ing on the actual metallicity. Planets are usually smaller in mass 
than brown dwarfs, but there is an overlapping region; a tiny por¬ 
tion of planetesimals, the so-called ‘super-planets’ may grow up 
to 20-40 Mjup by core-accretion (Mordasini et al. 2009), causing 
problems with separating these giant planets from brown dwarfs 
(Schneider et al. 2011). The maximum mass of brown dwarfs is 
about 75-80 Mjup (Baraffe et al. 2002). 

Chabrier et al. (2014) suggests that the borderline between 
giant planets and brown dwarfs should be linked to their different 
formation scenarios and not to the minimum mass required for 
deuterium ignition. Then the various formation and evolutionary 
mechanisms may be responsible for the nonequal frequency of 
brown dwarfs and planetary-mass objects around stars. 

While ~50% of solar-like dwarf stars have a stellar, ~30% 
have a low-mass (Earth- to Neptune-sized), and ~2.5% have 
a higher mass (Jupiter or bigger) planetary companion, only 
0.6-0.8% of solar-like stars have a brown dwarf within 5 AU 
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(Duquennoy & Mayor 1991; Vogt et al. 2002; Patel et al. 2007; 
Wittenmyer et al. 2009; Sahlmann et al. 2011; Dong & Zhu 
2013). The brown dwarf desert refers to this low occurrence rate 
of brown dwarfs as companion objects to main-sequence stars. 
However, the frequency of wide pairs consisting of a brown 
dwarf and a solar-like star (with a separation exceeding 5 AU) is 
significantly higher at over 2-3% (Ma & Ge 2014 and references 
therein). A more detailed overview of the different occurrence 
rates can be found in Ma & Ge (2014). 

Brown dwarfs may have a different formation mechanism 
from giant planets, which might form via core accretion (Alibert 
et al. 2005). Ma & Ge (2014) proposed that brown dwarfs below 
35 Mjup are formed from the protoplanetary disk via gravita¬ 
tional instability, while those with M > 55 Mjup form like stellar 
binaries via molecular cloud fragmentation. Between masses of 
35 and 55 Mjup there are a significant lack of brown dwarfs 
orbital periods shorter than 100 days (Ma & Ge 2014) as com¬ 
panions to stars. 

Alternative theories explaining the brown dwarf desert sug¬ 
gest that during the formation of a binary system consisting of 
a solar-type star and a brown dwarf, the migration process is 
so effective that brown dwarf companion spirals into the star 
and is engulfed. This would explain the paucity of very few 
brown dwarf companions found within 5 AU to solar-type stars 
(Armitage & Bonnell 2002). Combining this low frequency of 
close-in brown dwarfs with the geometric transit probability, 
transiting brown dwarfs should indeed be rare. 

In this paper we report the discovery of a transiting brown 
dwarf, CoRoT-33b. With a mass of 59 Jupiter mass, CoRoT-33b 
lies just below to the D-Li border for brown dwarfs at 65 Jupiter 
masses. We also give a preliminary estimate for the frequency of 
close-in brown dwarfs with orbital period less than ten days. 


2. Data 

2.1. CoRoT observations 

Flux measurements and transit detection 

CoRoT-33 (R = 14.25 magnitude) was observed by the CoRoT 
satellite (Auvergne et al. 2009; Baglin et al. 2007) in white light 
for 77.4 days between 8 July 2010 and 24 September 2010. In to¬ 
tal, 178,342 data points were collected. Various designations of 
the target are listed in Table 2 along with equatorial coordinates 
and magnitudes in different passbands. 

The first transits of CoRoT-33b were discovered in the so- 
called Alarm Mode which triggered an oversampling rate and 
spectroscopic follow-up observations. After 22 July 2010 the 
standard 512 s sampling rate was changed to 32 s (see Surace 
et al. 2008; Bonomo et al. 2012). In total, 13 transits of CoRoT- 
33b were observed. 

For the light-curve analysis we kept only those data points 
that were flagged with ‘0’ by the CoRoT automatic data-pipeline 
which is an indication of good measurement without comment. 
In total, we used 1626 data points obtained with the 512 s inte¬ 
gration time and 133,333 data points obtained with the 32 s in¬ 
tegration time for a total of 134,959 photometric measurements. 

The cleaned light curve normalized to its median value is 
shown in Fig. 1. Fortunately, no significant jumps, cosmic ray 
events, hot pixels, etc. affected this light curve. Also shown is a 
smoothed version of the data produced after applying a Savitzky- 
Golay filter. Variations due to rotational modulation can clearly 
be seen, which is studied in detail in Section [33] 
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Fig. 2. Finding chart and contamination source map for CoRoT- 
33. Red numbers denote stars whose contribution to the observed 
flux was taken into account; yellow numbers denote stars whose 
contamination was checked but was found negligible. Star with 
number 0 corresponds to CoRoT-33. 


Contamination 

Figure |7] shows the Palomar Observatory Sky Survey’s image 
of CoRoT-33 and its environment. The solid zigzag line repre¬ 
sents the CoRoT-photometric mask while the dashed line shows 
the boundaries of the CoRoT-imagette. CoRoT has a bi-prism in 
the optical pathway of its exoplanet channel so that a small, very 
low-resolution spectrum is obtained for brighter stars in the field. 
Therefore the point-spread function (PSF) of each object in the 
field of view of the exoplanet channel is roughly 46 x 23 arcsec- 
onds. We found that the contribution of stars outside the mask is 
very small, but there are two main contaminants inside the pho¬ 
tometric mask, denoted by Nos. 1 and 2. They are fainter than 
CoRoT-33 by 2.5 and 2.9 magnitudes. Because of the large PSF 
of CoRoT, the contaminating sources actually produce 13+4% 
of the total observed light in this mask. This contamination fac¬ 
tor was calculated using the procedure described in Pasternacki 
et al. (2011) and is consistent with the value found independently 
by Gardes et al. (2011). 


2.2. Radial velocity measurements 

The radial velocity (RV) foll ow-up of CoRoT-33 was started 
with the HARPS spectrograph (iMavor et al. II2003I) based on the 
3.6-m ESO telescope (La Silla, Chile) as part of the ESO large 
program 188.C-0779. To monitor the Moon background light on 
the second fibre, HARPS was used with the observing mode 
obj_AB, i.e., without simultaneous thorium-argon (ThAr). The 
exposure time varied between 25 minutes and 1 hour. A set of 
8 spectra was recorded for CoRoT-33 with HARPS between 28 
June 2013 and 15 August 2013. We reduced the HARPS data 
and computed RVs wit h the HARPS pipeline based on the cross - 
correlation technique (iBaranne et al. lll9^ iPepe et al. 1120021) . 
Radial velocities were obtained by weighted cross-correlation 
with a numerical K5 mask which yielded the smallest error bars. 
The signal-to-noise ratio (S/N) per pixel at 5500 A is in the range 
2 to 7. The first and the last measurements were affected and cor¬ 
rected from moonlight contamination. 

Pour additional RV measurements were acquired with the 
PIbre-fed Echelle Spec trograph (PIES; iPrandsen & Lindberd 
ll999l : lTilting et all2014l) mounted at the 2.56-m Nordic Optical 
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Table 1. Basic data of known transiting brown dwarfs, p is the mean density of the brown dwarf component. 


Name 

-^star/^O 

^star/^O 

[K] 

[Fe/H] 

P (days) 

e 

Mbd /Mjup 

Rbd/Rjup 

P [g/cm'] 

Ref. 

2M0535-05a'' 





9.779621(42) 

0.3225±0.0060 

56.7 ±4.8 

6.5±0.33 

0.26±0.06 

1 

2M0535-05b'' 





9.779621(42) 

0.3225±0.0060 

35.6 ±2.8 

5.0±0.25 

0.35±0.08 

1 

CoRoT-3b 

1.37 ±0.09 

1.56 ±0.09 

6740 ±140 

-0.02±0.06'' 

4.25680(5) 

0.0 

21.66 ± 1.0 

1.01 ±0.07 

26.4±5.6 

2 

CoRoT-15b 

1.32 ±0.12 

1 46‘^o.3i 

0 94®- 
008 

6350 ± 200 

±0.1±0.2 

3.06036(3) 

0 

63.3 ±4.1 

1 1 o-fO.30 

015 

59±29 

3 

CoRoT-33b 

0.86 ± 0.04 

5225 ± 80 

±0.44±0.10 

5.819143(18) 

0.0700 ± 0.0016 

59.0+1* 

1.10 ±0.53 

55±27 

4 

KELT-lb 

1.335 ±0.063 

1 4'7i-f0.045 

-P-P35 

1 ' 2 Q+ 0 .ri 
^•'^^- 0.10 

6516 ±49 

±0.052±0.079 

1.217513(15) 

0.007 

27.38 ± 0.93 

1.116t“* 

24.5ii, 

5 

Kepler-39b^ 

1 io+'’ 0 '' 

^■^^-0.07 

6260 ±140 

-0.29±0.10 

21.0874(2) 

'^•^^^-0.023 

18.00/«« 

20.1+11 

1 994-0.12 
^•^^- 0.10 

12.40/3J 

6 

Kepler-39b" 

1.40±0.10 

6350 ± 100 

±0.10±0.14 

21.087210(37) 

0.112 ±0.057 

1 944 - 0.09 
^•^^- 0.10 

13.0+^ 

7 

KOI-1896'' 

0.764 ±0.051 

0.733±0.017 

4952 ± 40 

-0.07±0.12 

30.3604467(5) 

0.2746 ± 0.0037 

78.0 ±3.4 

0.998 ± 0.023 

97.3±4.1 

8 

KOI-205b 

0.925 ±0.033 0.841 ±0.020 

5237 ± 60 

±0.14±0.12 

11.7201248(21) 

<0.031 

39.9 ± 1.0 

0.807 ± 0.022 

75.6±5.2 

9 

KOI-205b 

0 

0.04 

0.87 ± 0.020 

5400 ± 75 

±0.18±0.12 

11.720126(11) 

<0.015 

40.8=1’ 

0.82 ± 0.02 

90.9/1^ 

6 

KOI-415b 

0.94 ± 0.06 

L.LJ n ,0 

5810 ±80 

-0.24±0.11 

166.78805(22) 

0.698 ± 0.002 

62.14 ±2.69 

0 79+0-12 

157.4+^ 

10 

LHS 6343C" 

0.370 ± 0.009 0.378 ± 0.008 

3130 ±20 

±0.04±0.08 

12.71382(4) 

0.056 ± 0.032 

62.7 ± 2.4 

0.833 ±0.021 

109±8 

11 

WASP-30b 

1.166 ±0.026 1.295 ±0.019 

6201 ±97 

-0.08±0.10 

4.156736(13) 

0 

60.96 ± 0.89 

0.889 ±0.021 

107.6±1.1 

12 


Notes. References: 1; Stassun, Mathieu & Valenti (2006) ; 2M0535-05 2: Deleuil et al. (2008), 3: Bouchy et al. (2011a), 4: this study, 5: Siverd 
et al. (2012), 6: Bouchy et al. (2011b), 7: Bonomo et al. (2015), 8: Diaz et al. (2014), 9: Diaz et al. (2013), 10: Moutou et al. (2013), 11: Johnson 
et al. (2011), 12: Anderson et al. (2011) 

Notes, a: 2M0535-05 is an extreme young eclipsing system in which two brown dwarfs orbit each other, b: [M/H] value is reported in the 
reference. Notice that [M/77] « [Fe/7/]; we did not convert the inhomogeneous [Fe/H] to the same scale, c: aka K01-423b. d: Diaz et al. (2014) 
concluded that KOI-189b can be either a high-mass brown dwarf or a very low mass star, too, therefore its status is uncertain, e: the brown dwarf 
orbits companion A of a binary system, and data of the component A is given here. Star B has M = 0.30 ± O.OIMq, T^g = 3030 ± 30 K (Johnson 
etal. 2011) 


Table 2. IDs, coordinates, and magnitudes of CoRoT-33 from 
ExoDat (Deleuil et al. 2009) 


Designation 

CoRoT window ID 
CoRoTID 

2MASS 

USNO-A2 

USNO-Bl 

PPMXL 

Coordinates 

CoRoT-33 

LRc06_E2_1637 

105118236 

18383391±0537287 

0900-13338694 

0956-0378713 

5484010357803959995 


RA (J2000) 

18” 38“ 33.908'' 


Dec (J2000) 

±5° 37' 28.970" 


Magnitudes 

Filter 

Mag & Error 

Source 

B 

15.705±0.587 

mean of several 

sources 
(see ExoDat) 

R 

14.25 

USNO-A2.0 

1 

13.5 

PPMXL 

J 

13.238± 0.027 

2MASS 

H 

12.811± 0.026 

2MASS 

K 

12.707± 0.032 

2MASS 


Telescope (NOT) of Roque de los Muchachos Observatory (La 
Palma, Spain). The observations were carried out on 2, 3, 4, 
and 5 July 2014 under the CAT observing programme 79- 
NOT5/14A. We used the 1.3" med-res fibre, which provides 
a resolving power of R=s 47,000 in the spectral range 3600- 
7400 A. The weather was clear with seeing varying between 0.6 
and 1.4 " throughout the whole obs erving run. We followe d the 
observing stra t egy as described in iBuchhave et al.l (1201 Oh and 
iGandolfi et al.l (1201 5h and took three consecutive exposures of 
1200-1800 seconds per epoch observation to remove cosmic ray 
hits. We also acquired long-exposed (Texp = 22-24 s) ThAr spec¬ 
tra right before and after each epoch observation to trace the 


RV drift of the instrument. The data were reduced using stan¬ 
dard IRAF and IDL routines, which included bias subtraction, 
flat fielding, order tracing and extraction, and wavelength cali¬ 
bration. The S/N of the extracted spectra is 15-20 per pixel at 
5500 A. Radial velocity measurements were derived via S/N- 
weighted, multiorder, cross-correlation with the RV standard star 
HR 5777, and observed with the same instrument set-up as the 
target object. 

The HARPS and FIES RV measurements are listed in Table[3] 
along with their errors and the Julian date s of the observations in 
barycentric dynamical time (BJDtdb, see lEastman et alJl2010h . 


3. Analysis 

3. 1. Analysis of radial velocity data 

Full width half maximum (FWHM) and bisector span (BIS) 
of the cross-correlation function are also listed in Table 3 and 
show no significant variation in phase with the radial velocity. 
Furthermore the radial velocity peak-to-peak amplitude (~14 
km/s) is larger than the FWHM excluding a blended or back¬ 
ground binary scenario. 

An offset between the two instruments was calculated as 
part of the orbit fitting process and accounted for when com¬ 
bining the HARPS and FIES data sets. We first tried fitting a 
circular orbit to the RV data. Free parameters were the offset 
value for each data set, the true y-velocity and the amplitude (K) 
of the RV curve. The epoch and period were fixed to the val¬ 
ues of the photometric ephemeris (these were varied within the 
errors). This fit yielded a considerably large reduced ;if^-value: 
Amin “ allowed the epoch and period to vary in 

the fitting procedure. The resulting fit yielded a slightly smaller 
- 122. Interestingly, the period converged to a value that is 
compatible within Icr uncertainties of the period obtained from 
the photometric transits, but the resulting epoch obtained from 
the RV-analysis differed by 194 minutes from the photometric 
value. The error in the photometrically determined epoch is only 
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Fig. 1. Upper panel: full CoRoT light curve of CoRoT-33. The gray points represent the median-normalized raw data points after 
a 5-point width median filtering. We only used data points with flag “0” for this curve. The red line is a convolution of the raw 
light curve with a Savitzky-Golay filter that enhances the light-curve variations. Middle panel: Lomb-Scargle periodogram of the 
light curve of CoRoT-33. The horizontal dashed line denotes the 0.01% false-alarm probability (Scargle 1982). The vertical red line 
marks the rotation period of the star. Lower panel: autocorrelation function (ACF) of light curve, following the subtraction of the 
best fitting transit model. The red dashed line marks the peak corresponding to the rotation period of CoRoT-33 (see Section 3.5) 


6 minutes, so the epoch difference is significant at the 32cr level. 
Assuming a circular orbit the orbital fit to the RV data produced 
very large ^^-values or the RV data are incompatible with the 
photometric ephemeris. 

Therefore we tried fitting an eccentric orbit to the RV data. 
The results are quite satisfactory: ~ The orbit was 

found to have a slight eccentricity: e - 0.0700 + 0.0017. The 
results of this fit can be found in Table 4. 

Figure 3 shows the RV measurements - after correcting for 
the RV offset and subtracting the systemic radial velocity - phase 
folded to the orbital period (upper panel), along with the RV 
residuals (lower panel). 

Zakamska et al. (2011) found that a precise estimation of 
eccentricity requires that the radial velocity curve should have a 
S/N of 40, and the individual RV measurements should have a 
precision better than 1%. These requirements are fulfilled in our 


analysis. Their Eq. (7) gives an estimate of the precision of the 
eccentricity determination, i.e.. 


logcr(e) = 0.48 - 0.89 x log(/r Vv/cTobs) (1) 


where cr(e) is the expected precision on eccentricity, N - 12 is 
the number of the RV data points, CTobs = 0.046 km/s is their 
average (median) uncertainty. Substituting our values, we get 
cr(e) = 0.011 and thus our eccentricity determination and its 
error bar are reliable. The estimate of Zakamska et al. (2011) is 
the result of an average error estimate based only on RV, in our 
case the joint RV and light-curve fit yielded a smaller uncertainty 
range for the eccentricity (Sect. 3.4). The relatively large Kjcrohs, 
ratio enables us to achieve the aforementioned precision in the 
eccentricity. 
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Table 3. HARPS and FIES RV measurements of the CoRoT-33 system. FWHM is the average full-width at half maximum of the 
lines used (in Angstroms) and BIS is the corresponding bisector spans. S/N is the signal-to-noise ratio of the individual spectra. The 
FWHM of the FIFS CCFs are wider than those extracted from the HARPS spectra because of the different resolutions of the two 
spectrographs. 


Instrument 

BID 

Vrad [kms^‘ 

1 uncertainty [km s ' ] 

FWHM 

BIS 

S/N 

HARPS 

56471.70867 

27.159 

0.240 

10.8 

-0.08 

1.1 (moon) 

HARPS 

56473.65935 

14.793 

0.033 

10.1 

-0.04 

7.4 

HARPS 

56474.75692 

16.348 

0.127 

9.5 

-0.06 

2.0 

HARPS 

56475.74925 

23.725 

0.051 

9.9 

-0.07 

5.1 

HARPS 

56509.71203 

16.736 

0.058 

10.1 

0.03 

4.6 

HARPS 

56516.59952 

24.455 

0.040 

9.9 

-0.12 

6.0 

HARPS 

56518.67738 

25.197 

0.032 

9.9 

-0.06 

7.4 

HARPS 

56519.64904 

18.248 

0.085 

10.0 

0.04 

3.3 (moon) 

FIES 

56841.58071 

17.967 

0.026 

12.9 

0.014 

20 

FIES 

56842.62444 

25.419 

0.033 

12.7 

-0.056 

18 

FIES 

56843.60496 

28.120 

0.036 

12.5 

-0.044 

16 

FIES 

56844.61337 

24.856 

0.046 

12.6 

0.006 

15 
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Fig. 3. Upper panel; phase-folded RV measurements of CoRoT- 
33. Red circles represent the HARPS-measurements while green 
circles are the data points obtained by FIFS-instrument. Black 
solid line represents the eccentric orbit ht. The RV points and 
the ht is shifted by the y-velocity of the system. Lower panel: 
it shows the residuals of the ht. Vertical lines on the data points 
indicate their error bars. 


3.2. Spectral analysis of the host star 

We determine the spectroscopic parameters of the host star using 
the coadded HARPS and FIFS spectra. The S/N of the coadded 
data is relatively low, about 15 and 30 per pixel at 5500 A for 
HARPS and FIFS, respectively, owing to the low S/N of each 
individual spectrum 

We used the Spectroscopy Made Fasy package (SMF), 
version 4.1.2 (Valenti & Piskunov 1996, Valenti & Fischer 
2005) along with Atlas 12 or MARCS 2012 model atmospheres 
(Kurucz 2013, Meszaros et al. 2012) to determine the funda¬ 
mental photospheric parameters iteratively. SMF hts the ob¬ 
served spectrum directly to the synthesized model spectrum 
and minimizes the discrepancies using a nonlinear least-squares 
algorithm. The SMF utilizes input from the VALD database 
(iPiskunov et al.lll995l:lKuDka et al.l[T999l) . The uncertainties per 
measurement using SMF based on a sample of more than 1000 
stars was found by Valenti & Fischer (2005) to be 44 K for the ef¬ 


fective temperature Teff, 0.06 dex for the surface gravity logg», 
and 0.03dex for the metallicity [M/H]. This is in a statistical 
sense and somewhat optimistic. It was found from Valenti & 
Fischer (2005) that the scatter was systematically larger with an 
uncertainty in logg, of 0.1 dex and a scatter that occasionally 
could reach 0.3 dex. Fridlund et al. (2015, in preparation) have, 
in a systematic study of the CoRoT exoplanet host stars, con- 
hrmed these numbers. Furthermore, the errors in Teff are also af¬ 
fected by the problems with determining the shape of the Balmer 
lines accurately enough, as pointed out by Fuhrmann in a series 
of papers (see, e.g., Fuhrmann et al. 2011 and references therein; 
also Fridlund et al. 2015, in preparation). 

Torres et al. (2012) reported SME overestimates logg and 
correlates with T^ff and [Fe/H]; see also Brewer et al. (2015) 
that this issue is eliminated. Therefore, we carried out an in¬ 
dependent analysis of the spectrum using customized IDL soft¬ 
ware suite to ht the composite HAR PS and FIES spectra with a 
grid of theoretical m odel spectra from Castelli & Kuruczl (l2004l) . 
ICoelho et al.l (l2()()l . and iGustafsson et al.l (l2008l) . using spectral 
features that are sensitive to different photospheric parameters. 
Briehy, we used the wings of the Balmer lines to estimate the 
effective temperature of the star, and the Mgi 5167, 5173, and 


5184 A, the Cai 6162 and 6439 A, and the Na i D lines to de¬ 
termine logg*. The iron abundance [Fe/H] and microturbulent 
ve Mcity Vmkm was der i ved by applying the method described 
in [Blackwell & Shallia (Il979h on isolated Fe i and Fe ii lines. 
We adopted the calibra tion equations for Sun-like dwarf stars 
from lDovle et al.l (l2014t) to determine the macroturbulent veloc¬ 
ity, ''macro - The projected rotational velocity v sin i, was measured 
by htting the prohle of several clean and unblended metal lines. 

The two analyses provided consistent results well within the 
errors bars, regardless of the method and spectrum used. The h- 
nal adopted values - obtained as the weighted mean of the inde¬ 
pendent determinations - are listed in Table 4. We found Teff = 
5225 + 80K, logg, = 4.4 + 0.1 (cgs), vsini* = 5.7 + 0.4km/s and 
a c onsiderably high iron conten t: [Fe/H] = -1-0.44 + 0.10. Using 
the lStraizvs & Kurilienel(ll98lll calibration scale for dwarf stars, 
the effective temperature of CoRoT-33 translates to a G9 V spec¬ 
tral type. 


3.3. Estimation of stellar parameters 

The usual way to determine fundamental stellar parameters is to 
combine theoretical isochrones along with the measured mean 
density of the star, stellar metallicity and effective surface tem- 
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perature. The mean density can be obtained from transit dura¬ 
tion, the period and the impact parameter (e.g. Roberts 1899; 
Mochnacki 1981; Seager & Mallen-Ornelas 2003; Winn 2010). 
Stellar models are then selected that are able to reproduce the 
observed quantities. However, this method was found to be in¬ 
adequate in our case because we have a grazing transit (see 
Sect. 3.4), and as a consequence, a low-transit depth-to-noise 
ratio. 

Therefore the effective temperature, metallicity and log g, of 
the star were used to estimate the stellar mass and radius. We 
used the analytical stellar evolutionary tracks of Hurley et al. 
(2000) which have about 2% error relative to detailed numerical 
models, but are computationally fast. This yielded Mjtar = 0.86 + 
O.O4M0 and Rjtar = 0.94)^QQgRo for the stellar mass and radius, 
respectively. The age of the star is quite uncertain, but the object 
is definitely not young. This is supported by the fact that we did 
not find the presence of lithium in the spectrum. The lower age 
limit is 4.6 Gyrs and, most likely the star has an age of 11 Gyrs 
(see the derived stellar parameters in Table 4). 

As a sanity check we calculated the stellar radius indepen¬ 
dently of the stellar models. If the stellar rotational axis is per¬ 
pendicular to the orbital plane of the brown dwarf, then i, - 
85.5° (cf. Sect. 3.4). Assuming that the 8.936 days modulation 
of the light curve is the stellar rotational period (cf. Sect 3.5), 
then the stellar radius can be computed from 


Vrot sm 


27r7?5;tar sm 
Plot 


( 2 ) 


resulting in Rstar = 1.06 + 0.09Ro- This is well within the uncer¬ 
tainty range of isochrone-based value of Rstar = 0.94j^Q PgRo. 


3.4. Transit light-curve analysis 

To measure the transit parameters, we removed the stellar vari¬ 
ability from the light curve using the following procedure. First, 
we applied a median filtering to replace the outliers. We then fit 
the flux variations in a narrow vicinity of the out-of-transit data 
(separately from transit to transit) with a parabola. All points in 
and out of transit were divided by the corresponding parabola. 
The width of the window which defined the ‘vicinity of the tran¬ 
sit’ was found by visual inspection after assuming various trial 
lengths: ID, 1.5D, 2D, etc. up to 6D where D is the transit dura¬ 
tion. A +2D window around each transit center was found to be 
the most appropriate. 

The folded light curve using values binned by 201 s can be 
seen in Fig. 4. 

The modeling of the transit light curve is challenging be¬ 
cause of the faintness of the host star and small transit depth 
which result in a low S/N (transit depth/average noise level is 
~ 1). We used Transit Light Curve Modeling Code (TLCM, writ¬ 
ten by SzCs) which utilizes the Mandel & Agol (2002) model to 
describe the transit light-curve shape. A genetic algorithm was 
used to optimize the fit (Geem et al. 2001) and the error esti¬ 
mation was carried out using a simulated annealing chain of 10^ 
steps, starting from the best solution we found with the genetic 
algorithm procedure. 

When scaled semi-major axis (a/Rstar) is treated as a free pa¬ 
rameter, one gets the stellar mean density by rewriting Kepler’s 
3rd law, 

_ 37r I 

“ GP2(1 -H q) j 

where q - Mbd /iHstar is the mass ratio. We find the star would 
have pstar ~ 29 g cm“^ which is incompatible with the observed 


spectral type of G9V. The object G9V should have a mean den¬ 
sity of about 1.27 gcm“^. This fit also leads to 0.44 Jupiter-radii 
for the brown dwarf which results in a mean density of ~ 834 
gcm“^. This is unrealistically high as no brown dwarf is known 
to have a mean density higher 200 g cm“^ (cf. Table 1 and Ma & 
Ge 2014). These results are independent of how limb darkening 
is treated. In our case we used a quadratic limb darkening law 
fixed to values based on Sing’s (2010) tables. We also modeled 
the light curve allowing the limb darkening to vary as well as 
fixing one of the two limb darkening coefficients and allowing 
the other to be a free parameter. 

The transit depth is 0.28%, much shallower than the approx¬ 
imately 1% for a central transit of a Jupiter-sized brown dwarf 
around a late G-type star. The contamination by other stars in the 
photometric aperture does not explain the shallow transit depth. 
No additional contaminating source was revealed by the spec¬ 
trum of the host star. The transit curve is V-shaped and the tran¬ 
sit duration is just ~ 1.4 hours, much shorter than the expected 
~ 3 hours for a centrally transiting short-period substellar object. 
This suggests that we have a grazing eclipse and the transiting 
object crosses the apparent stellar disk on a chord shorter than 
the diameter of the star. 

The transit duration is (cf. Seager & Mallen-Ornelas 2003; 
Winn 2010): 


£)~ :^x —X V(1 +^)^' 

n a 


■ 


\ — e^ 


1 + e sin w 


( 4 ) 


where D is the measured transit duration (D = 1.4+0.1 hours), P 
is the orbital period, e is the eccentricity, w is the argument of the 
periastron, k is the radius ratio and b is the impact parameter. All 
parameters are known from photometry or from spectroscopy, 
except k and b. Since the radius ratio must be positive we have 
the requirement of k > 0. Substituting the measured values, we 
get that this requirement is fulfilled only if b > 0.91. This kind 
of a high impact parameter means we must be dealing with a 
grazing transit scenario. 

We conclude that the low S/N of the transit events prevents us 
from determining the transit shape with sufficient accuracy. We 
thus need to apply a penalty function (one may call a “prior”) to 
force the light-curve solution to converge to the right values: 


(p(fl/Ssiar)-1.269)2 

Q-X^+e- 2 x 0.2102 ( 5 ) 

and we minimized Q rather than of the light-curve fit. Here 
p represents the mean density of the star calculated for every 
iteration step of the fit via Eq. (3). The values 1.269 for the mean 
density and 0.210 for the width of the distribution are the mean 
density of the star and its uncertainty (in gcm“^), derived from 
the stellar mass and radius obtained in Section 3.2. 

In addition to the low S/N, the impact parameter is high, 
and it is difficult to fit the limb darkening coefficients if the im¬ 
pact parameter is larger than ~ 0.85. Since the brown dwarf 
only crosses a fraction of the stellar disk, the transit center is 
no longer flat-bottomed and the inner contact points disappear 
(Muller et al. 2013; Csizmadia et al. 2013). 

Another difficulty is connected to the fact that our theoretical 
knowledge of limb darkening is very poor close to the limb of 
the star. For instance, based on 38 Kepler light curves of tran¬ 
siting objects (12 of them are grazing), Muller et al. (2013) es¬ 
tablished that the quadratic limb darkening coefficients have a 
larger disagreement from the theoretically predicted values of 
Claret & Bloemen (2011) than the linear limb darkening coeffi- 
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cents 0 Difficulties in the theoretical limb darkening laws occur 
mostly at the limb because of the quadratic term. 

This is further illustrated by the fact that Claret & Bloemen 
(2011) and Sing (2010) used a plane-parallel stellar atmosphere 
model. When a more realistic spherically symmetric model is 
used, Neilson & Lester (2013) found a fast drop in the intensity 
close to the limb and the very edge of the star is predicted to be 
much darker than in Claret & Bloemen’s model (for comparison 
of the plane-parallel and spherically symmetric models, see Figs. 
2-5 of Neilson & Lester 2013). Since none of these new models 
have been checked against a large sample of eclipsing binaries 
and transiting objects covering a wide range of stellar tempera¬ 
tures, it is too early to accept these as a final description of the 
limb darkening law. In addition, these theoretical models do not 
take the stellar spots and faculae into account which can mod¬ 
ify the observable limb darkening coefficients (Csizmadia et al. 
2013). We should note that we can see evidence of the presence 
of spots on CoRoT-33, thus it is an active star. 

Therefore we carried out two light-curve fits with the Transit 
Light Curve Modeling Code (TLCM). During the fits, our free 
parameters were the scaled semimajor axis, the impact param- 
etefl the radius ratio and the epoch. The limb darkening com¬ 
binations u+ - Ua + Uh and u^ - Ua - Uh were fitted once (cf. 
Csizmadia et al. 2013 and Espinoza & Jordan 2015), and then 
for a second calculation they were fixed at theoretical values of 
Sing (2010). Eccentricity and argument of periastron could vary 
only within the limits of their uncertainties. 

In the case of the adjusted limb darkening coefficients, the 
values of the coefhcients are not constrained at all, but their un¬ 
certainties are large because of the low amplitude/noise ratio and 
the grazing transit. One cannot decide on a preferred solution 
based on the quality-parameter Q alone. Taking into account 
that the fixed limb darkening coefficient solution contains two 
less free parameters, we have chosen that solution for subsequent 
discussion in the paper. 

Then a joint fit of RV- and light-curve data by Exofast 
(Eastman et al. 2013) with adjusted limb darkening coefficients 
was carried out. This confirmed our solution in Table 4. The pa¬ 
rameters of this brown dwarf should be refined in the future us¬ 
ing better photometric measurements obtained with larger tele¬ 
scopes. 

The fit is shown in Table 4 and Eig. 4. 

3.5. Stellar variability analysis 

The light curve of CoRoT-33 exhibits periodic and quasiperi- 
odic flux variations with a peak-to-peak amplitude of about 3% 
(Eig. 1, upper panel). The Lomb-Scargle periodogram of the 
light curve is shown in Eig. 1. The strongest peak occurs at v 
= 0.1117 d^' (P - 8.95 d). Toward higher frequencies, there are 
two additional significant peaks at v = 0.2248 d“' (P - 4.44 d) 
and V = 0.3337 d ' (P = 3.00 d). The dominant peak is consistent 
with the expected rotational period, P^ot of the star calculated us¬ 
ing the spectroscopically measured rotational velocity and stel¬ 
lar radius. Therefore, this signal is most likely due to spots since 
CoRoT-33 should have a modest level of magnetic activity given 


^ Their limb darkening formula that we also use is / = -g)- 

U 2 ( l - p)^, where Ig is the intensity at the center of the apparent stellar 
disk, Ml, U 2 are the linear and quadratic limb darkening coefficients, 
respectively; p = cos y; and y is the angle between the line of sight and 
the surface normal vector at the stellar surface point. 

^ We used b = t-^-4 — to take the effect of eccentric orbit on the 

Jfsiar 1+esmw 

impact parameter into account. 


Table 4. Physical and geometrical parameters of the CoRoT-33 
system. Inclination (i) was calculated from the a/Rstar ratio and 
from the impact parameter i>. The parameter can be cal¬ 

culated from the orbital period and from the afRstar value (see 
e.g. Winn 2010). 


Determined from photometry 

Epoch of transit Tg [BJD-2450000] 

6676.3992±0.0037 

Epoch of periastron Tg [BJD-2450000] 

6677.7130±0.0140 

Orbital period (days) 

5.819143±0.000018 

Duration of the transit (hours) 

1.4 hours 

Depth of the transit (%) 

0.28% 

Determined from RV measurements 

Orbital eccentricity e 

0.0700 ± 0.0016 

Argument of periastron cj [deg] 

179.3 ±0.87 

RV semi-amplitude K [kms“‘] 

7.0609 ± 0.0094 

Systemic velocity Vy [kms^*], HARPS 

21.5375 ±0.0158 

Systemic velocity Vy [kms^*], FIES 

21.5544 ± 0.0222 

O-C residuals" [ m s“*] 

36 

Determined from spectral analysis of the star 

Teff [K] 

5225±80 

logg. [cgs] 

4.4±0.1 

[Fe/H] 

0.44±0.1 

[Ni/H] 

0.4±0.1 

[V/H] 

0.4±0.1 

[Mg/H] 

0.4±0.1 

[Ca/H] 

0.3±0.1 

[Si/H] 

0.2±0.1 

vsini, [kms^'] 

5.7±0.4 

Spectral type 

G9V 

V ■ 

* mic 

0.86±0.1 km/s 

V’’ 

mac 

2.7±0.6 km/s 

Determined from light curve modeling 

^/^star 

13.23±1.17 

b 

1.04±0.06 

^planet [dcg] 

85.5±0.5 

k 

0.12±0.04 

contamination [%] 

13±4" 

Combined results 

Stellar mass [solar] 

0.86!!;;- 

Stellar radius Rst [solar] 

0-94!°:'‘^ 

Stellar age [Gyrs] 

> 4.6 

Orbital semi-major axis a [AU] 

0.0579“' 

Brown dwarf mass Mbd [Mj ] 

59.0+1* 

Brown dwarf radius Rbd [R j ] 

1.10 ±0.53 

Brown dwarf mean density pbd [g cm“^] 

55 ±29 


Notes. FI Root mean square of the residuals of the RV curve. F) Fixed at 
this value, based on the calibration by Doyle et al. (2014). F) Calculated 
from period and masses via Kepler’s third law, not from RV. From 
modeling results. 


its relatively fast rotation rate. Interestingly, this rotational period 
is ~ 3/2 X Porb- The other two peaks correspond to the first and 
second rotational harmonics (Piot/2, Riot/3, respectively). This 
probably reflects a complex pattern on the stellar surface. Eor in¬ 
stance, Prot/2 can result from two spot groups on opposite sides 

of the star. _ 

Eo llowing the g uidelines describ e d in iMcOuillan et al.l 
(12013h and revised in iMcOuillan et alJ (l2014l) . we use the au¬ 
tocorrelation function (ACE) to confirm these results. The ACE 
of the CoRoT-33 light curve shows a strong correlation peak at 
the rotational period of approximately nine days, consistent with 
the Lomb-Scargle periodogram, followed by a sequence of addi¬ 
tional local maxima at integer multiples of this value. This is the 
result of the viewing geometry which repeats after multiple rota- 


7 
























Csizmadia et al.; CoRoT-33b: a synchronized brown dwarf 


Table 5. Results of three light-curve solutions of CoRoT-33. Idc stands for limb darkening coefficients. Inclination is calculated 
from the scaled semi-major axis, impact parameter, eccentricity and argument of periastron. 


Parameter 

fixed Idc (TLCM) 

adjusted Idc (TLCM) 

ExoFast 

^/^star 

13.23±1.17 

13.24±0.33 

14.20±1.25 

b 

1.04±0.06 

1.00±0.05 

1.053±0.065 

^planet [^^§1 

85.5±0.5 

85.7±0.4 

85.70±0.53 

k 

0.12±0.04 

0.078±0.034 

0.139±0.039 


0.57 

0.13±0.36 

- 

U- 

0.12 

0.33±0.99 

- 

Ml 



0.538±0.053 

U2 



0.171±0.051 

contamination [%] 

13±4 

13±4 

13% (fixed) 

Q 

1.257 

1.257 




- 0.02 - 0.01 0.00 0.01 0.02 
Phase 


Fig. 4. Light-curve solution (solid red line) and phase-folded, 
binned data points (black circles) of CoRoT-33b transits. Lower 
panel shows the residuals of the fit. Note that transit depth is only 
0.28% on an R = 14.25 magnitude star. 


tions (Fig. 1, lower panel). We estimate the position of the peaks 
by fitting Gaussian functions to the sequence of ACF maxima, 
and define the rotation period as the slope of a straight-line fit 
to the ACF peak positions as a function of peak numbers. This 
yields a stellar rotation period of Prot = 8.936 + 0.015 days. 

We also investigate the periodicities in the light curve using 
the package MUFRAN developed by Kollath (1990). This soft¬ 
ware is an efficient tool for detecting periodic patterns in time 
series and is based on the Fourier transform. 

Owing to temporal gaps in the data some false frequencies 
(aliases) can appear in the Fourier power spectrum. These alias 
frequencies are centered on the real signals offset from those 
peaks as indicated by the spectral window function of the Fourier 
transform. The spectral window (Fig. |5]l indicates a negligible 
aliasing caused by the data sampling. 

The Fourier amplitude spectrum of the whole data set is 
shown in Fig. |6] for this wide interval. The multiple peaks near 
14 c/d are artifacts due to scattered earthshine arising from the 
103 min orbit of the CoRoT satellite. For better visibility, the 
low-frequency part of the amplitude spectrum is shown in the 
insert of Fig. |6] This Fourier analysis confirmed the 8.936 day 
rotational period and the presence of its harmonics. The period 
analysis was carried out separately for the first and second half 


of the data set, too, and it showed that the main frequencies do 
not change with time. 

3.6. Search for beaming effect 

One can ask whether one of the frequencies in the Fourier spec¬ 
tra is related to beaming effect (Zucker et al. 2007). Mazeh & 
Faigler (2010) cleaned and binned the light curve to 100 points 
to see this effect in CoRoT-3. Unfortunately, we cannot follow 
their approach. We can estimate the magnitude of this effect us¬ 
ing their formula: 

F(t) - F(i+A^sin(Lj(t-to))-A^os(io(t-tQ))-AeCos{2tj(t-to))(6) 

where w = 2;7r/RorbitaU and Ab, Ai , Ag are the amplitudes of the 
beaming effect, the variation coming from reflected light and 
from the ellipsoidal shape of the star. All these effects act on 
the same time-scale, therefore we see a periodic signal at the 
orbital period and its half value. Using the formulae of Mazeh 
& Faigler (2010), we expect Ab = 94 ppm, Ar = 9 ppm, and 
Ae = 54 ppm for a total peak-to-peak amplitude of 231 ppm. 
Since the eccentricity is small, we can neglect its effect. Notice 
that these estimates are uncertain by several percent (Mazeh & 
Faigler 2010). 

For an albedo of 0.1 and 0.9, the equilibrium temperature of 
the brown dwarf object would be between 1400-800 K in this 
system. Baraffe et al. (2003) predicts a surface temperature of 
1100 - 1200 K for a brown dwarf with a mass and age consistent 
with our object. Thus the brown dwarf contributes only 0.2 ppm 
to the observed visible light so it can be neglected. 

Using the approximate formula of Aigrain et al. (2009), we 
get that the two hours average noise level of CoRoT at the mag¬ 
nitude of CoRoT-33 is about 365 ppm (not taking the aging of 
the CCD detector into account). This means that we have to ap¬ 
ply very strong binning to see the effect. To detect this signal 
at a S/N of 3 we had to apply a binning of two days. Therefore 
we choose a different approach from that of Mazeh & Faigler 
(2010); namely we utilize the Fourier-spectrum. 

The cleaned light curve was taken again and a five-point 
median-filtering was carried out. We discovered a significant 
peak at 5.93 + 0.32 days with an amplitude of 318 ppm in the 
Fourier-spectrum after a cleaning process. This peak is signif¬ 
icant at the 13cr level (its significance was estimated by using 
Eq. (21) of Kjeldsen & Frandsen 1992). Its amplitude is 50% 
bigger than our expectation for the peak at the orbital period, but 
the expectation is also uncertain by several percent. Although 
this peak is quite close to the orbital period, it can also be identi¬ 
fied as a harmonic of the rotational period of the star; the period 
ratio of this peak and the stellar rotational period is exactly 1.5. 


0.002 - 

0.001 - 

0.000 - 

- 0.001 - 

- 0.002 - 
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Fig. 5. Window function of CoRoT-33. 
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Fig. 6. Fourier amplitude spectrum of the light curve of CoRoT- 
33. Frequencies are in cycle/day units. The insert shows a zoom 
into the low-frequency part. 


Further study is required to separate the stellar rotational modu¬ 
lation and the beaming effects. 

3.7. A search for occultations 

A brown dwarf in a close-in orbit may be occulted by the parent 
star (e.g. Winn 2010) which could be helpful to constrain the im¬ 
pact parameter and other properties of the system with higher ac¬ 
curacy. We searched for these kinds of occultations of the brown 
dwarf in the light curve of CoRoT-33 using a Bayesian model 
selection (e.g. Kass & Raftery 1995, Gregory 2010), wich are 
similar to the approaches used in Parviainen et al. (2013, 2014) 
and Gandolfi et al. (2015). The stellar variability was taken into 
account using Gaussian processes (Roberts et al. 2013; Gibson 
et al. 2011; Rasmussen et al. 2006). The nonzero orbital eccen¬ 


tricity means that the impact parameter for the occultation, bo, 
can be significantly different than the impact parameter for the 
transit, bt, or there might be no occultation. The technical details 
of this search are given in the Appendix. Our analysis was unable 
to find the occultation signal. The missing occultation signal can 
also be consistent with Baraffe et al.’s (2003) models. 

4. Tides and stellar rotational properties 

The rotation period is small for a G9V star of its age. The brak¬ 
ing associated with the activity of a G9V star is not as efficient as 
it should be if we take the single-star scenario of Bouvier et al. 
(1997) or Mamajek & Hillenbrand (2008). Because of the heavy 
companion orbiting the star at a short distance, the tidal torques 
on the star are strong enough to accelerate the stellar rotation. 
Several factors concur to make this system one of the best suit¬ 
able for the study of the interplay of magnetic braking and tidal 
evolution; the high mass of CoRoT-33b, the short distance from 
it to the star, and the age of the system. 

The simulations done by Ferraz-Mello et al. (2015) show 
that in a system like CoRoT-33, the magnetic braking may be 
very efficient in the beginning and this drives the rotational pe¬ 
riod to a value somewhat larger than the current observed stellar 
rotational period, when braking and tidal evolution equilibrate 
themselves. Subsequently, the system evolves losing energy, and 
the companion orbit spirals down toward the star and both the or¬ 
bital period of the companion and the rotation period of the star 
slowly decrease for the remaining life of the system, and it can 
possibly reach such commensurability just by chance. The fact 
that the orbital period and the stellar rotation period are strik¬ 
ingly close to the 2/3 ratio is not predicted in usual spin-orbit 
dynamical theories. Beky et al. (2014) listed six systems with 
a hot Jupiter where the planet and the star exhibit similar syn¬ 
chronization between these periods. They consider that the stel¬ 
lar differential rotation profile may happen to include a period 
at some latitude that is commensurable to the planetary orbit. 
These influence the displacement of the star features responsible 
for the periodic variation of the light of the star. The quality of 
our characterization of CoRoT-33 is in favors of it playing a key 
role in the study of these sorts of interactions and their influence 
on the tidal evolution of the system. 

5. Conclusions 

We report the detection of a transiting 59.0/;[ * Jupiter-mass 
brown dwarf that orbits a G9V star in ~5.82 days. This is the 
10th transiting object in the so-called brown dwarf desert not 
counting the transiting double brown dwarf system 2M0535- 
05 here). This desert is starting to be populated by more and 
more objects in recent years. The radius of the brown dwarf is 
1.10 + 0.53 Jupiter radii and its mean density is 55 + 29 gcm“^. 
The object is close to the 65 Mjup limit which separates the 
lower mass, deuterium-burning brown dwarfs from the higher 
mass, lithium-burning brown dwarfs. The host star seems to be 
an evolved, old, metal-rich star. The radius of the brown dwarf 
is not known to high precision because of a grazing transit. On 
the other hand, the mass is measured with higher accuracy from 
the RV observations, and its mass uncertainty is dominated by 
the uncertainty in the stellar mass. Better photometric precision 
using larger telescopes is required to improve the radius and thus 
density determination. 

The statistical analysis of 65 brown dwarf companions to 
stars with orbital period less than 100 days by Ma & Ge (2014) 
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a [AU] 


Fig. 7. Occurrence rate of brown dwarfs as companions to solar¬ 
like stars vs orbital separation (a). Horizontal lines denote the 
found occurrence rates and the orbital separation range investi¬ 
gated while vertical lines denote the error bars. Different colors 
mark different studies; black is Metchev & Hillenbrand (2009); 
blue is Lafreniere et al. (2007), Grether & Lineweaver (2006), 
orange is Sahlmann et al. (2011), brown is Wittenmyer et al. 
(2009), violet is Patel et al. (2007), red is this study. These 
rather rough estimates indicate a tendency; the occurrence rate 
increases with increasing orbital separations. 

showed that their host stars are not metal rich stars. This may be 
related to a different formation scenario of brown dwarfs com¬ 
pared to giant planets. The host stars of giant planets seem to 
be more metal rich (Johnson et al. 2010). However, CoRoT- 
33 is metal rich with [Fe/H] = -1-0.44 + 0.10 and its abun¬ 
dance may exceed the previous record holder, HAT-P-13 with 
[Fe/H] = 0.41 + 0.08 (Bakos et al. 2009; Ma & Ge 20140. 

The light curve of CoRoT-33 shows a rotational modulation 
with a period of 8.936 d with a peak-to-peak amplitude of ~3%. 
This measured rotational period agrees with the calculated value 
based on the vsini, measurement and its stellar radius from 
isochrones (assuming i, « ibd)- The Fourier spectrum is com¬ 
plex, making the target ideal for future studies of spot activity 
on this star. In particular, multicolor photometry using large tele¬ 
scopes can help to determine the spot temperatures. 

The rotational period of the star is very close to the 2:3 com- 
mensurability with the orbital period of the companion brown 
dwarf and this system is an interesting test case for checking and 
calibrating theories of the physical interactions between one star 
and a close-in companion. This relatively old system shows an 
eccentric orbit; the study of the tidal evolution of this system 
(Ferraz-Mello et al. 2015) shows that because of the distance 
between the brown dwarf and the star (0.0626 AU), an existing 
initial eccentricity is not damped to zero during the stellar life¬ 
time, however, the presently observed rotational period of the 
star cannot be explained without taking the interplay between 
magnetic braking of the star and tidal forces into account. 

One can assume that close-in hot Jupiters and brown dwarfs 
have the same transit detection bias since they both have ap¬ 
proximately the same radius. Therefore, the ratio of the num¬ 
ber of these detected systems should reflect the respective num¬ 
ber frequencies of these objects. Until now, three brown dwarfs 
(CoRoT-3b, -15b and 33b) have been detected from the CoRoT 
data, and 21 hot and normal Jupiters (CoRoT-lb, -2b, -4b, -5b, 

The brown dwarf HAT-P-13c is not transiting, it was detected by 
RV (Bakos et al. 2009). 


-6b, -9 —14b, 16b —21b, -23b, 25b —29b) most with peri¬ 
ods P < 10 daysQ The relative frequency, based on the CoRoT- 
sample and hence the observational biases are removed, is of a; 
14%, but one has to take into account that so far we have only a 
small sample statistic. 

If we take all transit surveys, then the brown dwarf/hot 
Jupiter ratio falls down to 0.05%. For instance, the WASP survey 
has detected one brown dwarf (see Table 1) and 96 hot Jupiters. 
HAT-P-survey has detected 33 hot Jupiters and none with a mass 
exceeding 8 Mjup; there is no transiting brown dwarf in the 
HAT-P-sample (see also footnote 4). It is not well known how 
the targets for follow-up and the stellar samples are selected in 
these surveys. Better statistics should come from Kepler once the 
follow-up measurements for the transit candidates are completed 
which may resolve the discrepancy between ground-based and 
CoRoT observations. A major factor can be that ground-based 
surveys are biased for the detection of transits smaller than 1%. 

The true occurrence rate of hot Jupiters is 1.2 + 0.4% around 
FGK dwarf stars for periods P < 10 days (Wright et al. 2012). 
Thus, the relative frequency of brown dwarfs to hot Jupiters 
can be scaled to the true frequency which means that the ac¬ 
tual occurrence rate of brown dwarfs as companions to FGK 
dwarfs would be ~ 0.2% in the short period (P < 10 days) 
range. This is much smaller than the frequency rate of brown 
dwarfs as FGK companions at larger distance (see Introduction). 
One can plot the different brown dwarf occurrence rate esti¬ 
mates mentioned in the Introduction and determined here as a 
function of the star-companion distance. The result is shown in 
Fig. 7. A fit to those data yielded that the brown dwarf occur¬ 
rence rate (/) around FGK dwarfs can be characterized roughly 
as / = 0.55/;Qgj(a/lAU)°'^^=^‘’°®, where a is the semimajor axis 
of the orbit. This result is not robust yet because the occurrence 
rate estimates suffer from small number statistics, but a tendency 
might be visible already. More observational studies are needed 
to establish the occurrence rate-orbital separation relationship 
for solar-like star+brown dwarf pairs. 

The early indication is that the occurrence rate of close-in 
brown dwarf (P < 10 days) is six times smaller than that of 
hot Jupiters in the same period range. It is not clear whether 
this relative occurrence rate is a consequence of the primordial 
conditions (i.e., much fewer brown dwarfs were formed in the 
protoplanetary disk than giant planets), or is caused by the higher 
efficiency of engulfing the companion by the host stars due to 
fast spiralization, as Armitage & Bonnell (2002) presumed. It is 
also possible that observational biases act: for instance, CoRoT’s 
stellar sample or its follow-up strategy or its higher photometric 
precision than ground-based surveys were simply more sensitive 
to brown dwarfs. 

For larger sized planets, the planet-frequency decreases for 
shorter periods (Dong & Zhu 2013). The brown dwarfs seems to 
follow the same pattern, too, but we think there is an indication 
that the frequency of close-in brown dwarfs drops more steeply 
than that of close-in giant planets. However, further detection 
studies are needed to establish the true frequency of close-in 
brown dwarfs. 
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Appendix A: Bayesian modei seiection 


Here we describe our methodology which was used for search¬ 
ing the possible occultation of the brown dwarf. 

We consider two models, Mq without an occultation signal 
and Ml with an occultation signal. We also use the Bayes factor 
Bio, the ratio of the Bayesian evidences Zm, to assess whether 
the occultation model is significantly preferred over the non- 
occultation model. The evidence Zm for a model M is calcu¬ 
lated by marginalizing the posterior probability over the model 
parameter space. 


Zm - 


J P(fflD)d0^ J 


P(6l|T>)d6»= F(0)F(DI0)d0, 


(A.l) 


where 0 is the model-specific parameter vector, D is the observa¬ 
tional data, F(0ID) is the posterior probability density, F(0) is the 
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prior density, and P{D\0) is the likelihood for the observational 
data. 

We assume the photometric noise to be normally distributed, 
and express the likelihood as 

P{D\e) = - ^ («D In In + In det S + , (A.2) 

where no is the number of data points, S is the covariance ma¬ 
trix, and r is the observed-measured residual vector. The covari¬ 
ance matrix elements are dehned by the Gaussian process (GP) 
covariance kernel, k{xi, xj, 0), where x are input parameter vec¬ 
tors for each data point, and 0 is a kernel hyperparameter vector. 

We only use one GP input parameter, mid-exposure time. We 
decided not to marginalize over the GP hyperparameters (we hx 
them to values optimized to the light-curve), but we did decide 
to use two different GP kernels to assess the sensitivity of our 
analysis on the choice of kernel. The two kernels used were the 
squared exponential (SE) kernel 


ksE(ti, tj, h, d) = exp 




(A.3) 


180 ° 



and the exponential (E) kernel 

tj, h, d) = exp 

where t is the mid-exposure time, h is the output scale, and d 
is the input scale. The SE kernel leads to inhnitely-differentiable 
smooth functions, while the E kernel leads to once-differentiable 
functions allowing for sharper changes (which is a more realistic 
choice considering the noise properties of CoRoT light curves). 

Parameterization and priors 

We construct the priors for the orbital parameters and radius ratio 
based on the corresponding marginal posteriors from the transit 
and RV modeling, which leaves the surface flux ratio as the only 
truly unconstrained parameter. 

In parallel with the main model selection analysis, we carry 
out a more explorative occultation search by mapping the Bu)- 
space as a function of sliding prior on w. This corresponds 
roughly to making n model comparisons for a set of proposi¬ 
tions, each differing in the prior set on oj. We use an uniform 
prior dehned by its center and width, and let the prior center 
slide from 0 to 2n. 


Fig.B.l. Differences between the Mq and Mi log posterior sam¬ 
ples with the SE kernel (corresponding to the log posterior ra¬ 
tios) for a set of 15000 posterior samples mapped as a function 
of argument of periastron (angle). A uniform prior from 0 to 2n 
has been set on w, and all the parameters have been drawn from 
their corresponding priors. The shown value, sample log pos¬ 
terior difference, does not correspond to the Bayesian evidence 
Z, but is used as an explorative tool. The radial spread in the 
difference between to values of n and 2n is explained by the or¬ 
bital geometry. The occultation signals are stronger within this 
interval than on the other half, and the data is able to exclude 
these signals leading to small log posterior ratio. The occulta¬ 
tion signals between cu of 0 and n are too faint to be excluded by 
the data, and the with-occultation corresponds to the occultation 
model within uncertainties. 


(A.4) 


Appendix B: Results 

We did not hnd support for the occultation model over the no- 
occultation model when we constrain to using a prior based on 
RV and transit modelling. The geometry, to ~ 180°, combined 
with the high impact parameter makes the detection of even a 
very strong signal unlikely. 

The mapping of Z?io as a function of to results in a tentative 
occultation signal candidate near 270°, as shown in Eig. lB.ll The 
result is curious, since this is the only geometry where an occul¬ 
tation could in theory be detected, while a false signal of instru¬ 
mental or astrophysical origin could present itself anywhere in 
w-space. However, the small number of orbits covered combined 
with the faintness of the signal candidate makes this impossible 
to confirm. If the occultation signal were to be real, it would 
be discrepant with the RV observations. This could be allevi¬ 
ated if the RV observations were to contain an unknown noise 


source, but the increase in the uncertainty per RV observation 
would need to be around 0.2-0.3 km/s. 
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The two GP kernels yield very similar results, meaning our 
analysis is not sensitive to the choice of the GP kernel. 
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